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ABSTRACT 
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Although measuring the gas metallicity in galaxies at various redshifts is crucial to constrain galaxy evolutionary sce- 
narios, only rest-frame optical emission lines have been generally used to measure the metallicity. This has prevented 
us to accurately measure the metallicity of dust-obscured galaxies, and accordingly to understand the chemical evolu- 
tion of dusty populations, such as ultraluminous infrared galaxies. Here we propose diagnostics of the gas metallicity 
based on infrared fine structure emission lines, which are nearly unaffected by dust extinction even the most obscured 
systems. Specifically, we focus on fine-structure lines arising mostly from Hll regions, not in photo-dissociation regions, 
to minimize the dependence and uncertainties of the metallicity diagnostics from various physical parameters. Based 
(— I ' on photoionization models, we show that the emission-line flux ratio of ([Oni]51.80-|-[Oin]88.33)/[Nni]57.21 is an ex- 

Q^l cellent tracer of the gas metallicity. The individual line ratios [Oni]51.80/[Nni]57.21 or [Olll]88.33/[Nin]57.21 can also 

I . be used as diagnostics of the metallicity, but they suffer a stronger dependence on the gas density. The line ratios 

O ' [Oni]88.33/[Oni]51.80 and [Nn]121.7/[Nni]57.21 can be used to measure and, therefore, account for the dependences 

^ ' on the of the gas density and ionization parameter, respectively. We show that these diagnostic fine-structure lines are 

, detectable with Herschel in luminous infrared galaxies out z ~ 0.4. Metallicity measurements with these fine-structure 

, lines will be feasible at relatively high redshift (z ~ 1 or more) with SPICA, the future infrared space observatory. 

Key words, galaxies: abundances - galaxies: evolution - galaxies: general - galaxies: ISM - Hll regions - infrared: 
galaxies - infrared: ISM 
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1. Introduction such rest-frame optical diagnostics cannot probe the metal- 

*~v^i . ... „ , licity of the regions affected by significant dust extinction, 

N The metallicity of gas and stars m galaxies is one of the ^^-^^ ^^^^ ^^^^^ ^^^^ ^^^^i^g galaxies, especially 

(N most important properties to distinguish various galaxy ^-^-^ redshift. This is because the fraction of dusty galax- 

evolutionary scenarios, since metals result from the cumu- ^^^^ ultra-luminous infrared galaxies (ULIRGs) in- 

O lative star-formation activity and gas inflow/outflow his- ^^^^^^ ^ f^^^^i^^ of redshift (e.g., Le Floc'h et al. 2005). 

7^. ; tory m galaxies Since the gas metallicity of ga axies is m ^^^^^^ ^^^-^^^ metallicity studies probably probe only the 

> , most cases much easier to measure than the stellar metal- ^^^^^^ j^^g extinguished region of star forming galaxies, 

- - hcity, past studies on the metallicity of galaxies and its ^^-^^ ^^^-^^^ substantially from the metallicity of in- 

eyolution have mostly focused on the metallicity of the gas ^^^^ ^^^-^^ ^^^^ obscured regions of star forma- 

phase (e.g., McCah et al. 1985; Zaritsky et al 1994; Pettim ^-^^ ^ ^j^^^ indication of such a mismatch was recently 

et al. 2001). Several extensive and detailed studies have ^^^^ submilhmeter galaxies by Santini et al. 

been performed on the gas metallicity m local galaxies, the (2010) , who used dust-mass measurements to infer an met- 

connection with other galaxy properties and the metallic- ^^^^^^^ ^^^^^ magnitude higher than inferred from 

oT/*o°'' !^onn°.'''^i'u''P°''^' tt ""r^' ''^ Optical metallicity diagnostics. If such mismatch applies to 

: onnf ' ° , °onAn nn?n^'' o ° a significant fractiou of actively star forming galaxies, this 

al. 2008; Mannucci et al. 2009, 2010). Such observational ^^^j^ -^^^^ ^ ^^-^^ ^^^-^-^^ p^^^ ^^^^l^^ 

studies have provided strong constraints on evolutionary ^metallicity evolution in galaxies based on optical diagnos- 

scenarios (e g., de Rossi et al 2007; Kobayashi et al. 2007; ^-^^^ ^^thin this context we note that, even for starbursts 

Brooks et al. 2007; Fmlator & Dave 2008). ^l^i^h have been well-studied in optical, a combination of 

The vast majority of previous studies on the gas metal- ^^^^^ infrared measurements has conclusively demon- 

hcity m galaxies used rest-frame optical diagnostics (see, gyrated that some of the most active star-forming sites are 

e.g., Nagao et al. 2006a and references therem). However, optically obscured (e.g., Gorjian et al. 2001; Vacca et al. 
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2002; Soifer et al. 2008). 
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Fig. 1. SEDs of the ionizing continuum adopted in our pho- 
toionization models, which are starburst99 constant star- 
formation model spectra (Leitherer et al. 1999) with an 
age of 1 Myr and an a = 2.35 IMF (Miow = lAf© and 
Map = IOOMq). Although most of our analysis adopt these 
spectra, models with different ages (3, 10, and 30 Myr) are 
investigated in §3.3. The blue, cyan, green, orange, and red 
spectra denote a model with Zgtar = 0.05Zq, 0.2Zq, OAZq, 
1.0^0, and 2.0Zq, respectively. The Zstar = 2.0Zq spec- 
trum is used also for our Zgas — 3.0Zq calculations (see 
the main text). All spectra are normalized at 5000A (~ 
0.18 Ryd). 



The access to far-infrared diagnostics of the gas metal- 
licity would overcome the dust extinction problems plagu- 
ing optical measurements. Observational studies on gas 
metallicity exploiting infrared indicators have already been 
carried out by using Infrared Space Observatory (ISO). 
Verma et al. (2003) investigated mid-infrared spectra of 12 
nearby starburst galaxies and found a strong correlation 
between gas metallicity and gas excitation, which is also 
seen in some optical studies (e.g., Nagao et al. 2006; Maier 
et al. 2006). Garnett et al. (2004) studied chemical proper- 
ties of gas in M51 based on the ISO data and found that 
the C/0 abundance ratio in M51 is consistent with the 
solar neighborhood value. This infrared-based metallicity 
studies is expected to progress dramatically by using the 
on-going and forthcoming observational facilities such as 
Herschel, JWST, SPICA, and ALMA. Within this context, 
it is important to investigate in detail model predictions 
yielding calibrations between the gas metallicity and flux 
ratios of emission lines at long wavelengths. Motivated by 
this, we present the results of photoionization model cal- 
culations and the theoretical calibrations of gas metallicity 
diagnostics based on fine-structure emission lines in mid- 
and far-infrared wavelength ranges. 



2. The method 

In the mid-infrared to sub-millimeter wavelength range, 
there are many fine-structure emission lines radiated by 
various ions in various excitation levels. Although some of 
them arise within Hii regions, some others arise in photo- 
dissociation regions (PDRs) that surrounds the inner Hii 
regions. For instance the [Cii] 157.6 emission, that is the 
strongest emission line in spectra of most galaxies and thus 
detected even in very distant galaxies (e.g. Maiolino et 
al. 2005; lono et al. 2006; Maiolino et al. 2009, Hailey- 
Dunsheath et al. 2010), is mainly from PDRs rather than 
Hii regions (e.g., Crawford et al. 1985; Stacey et al. 1991). 
Some emission lines arise both in Hii regions and PDRs 
(e.g., Petuchowski & Bennett 1993; Heiles 1994; Aannestad 
& Emery 2003; Abel 2006). The latter class of hues are the 
most difficult to model, since they have a more complex de- 
pendence on the various physical parameters of the clouds, 
and therefore they are not ideal metallicity diagnostics. 
Even lines arising mostly from PDR are not optimal to be 
used as metallicity diagnostics, since modeling PDRs is very 
challenging, because involves complex physical mechanisms 
of the gas heating and of radiative transport. The modeling 
of Hii regions is much more easy and generally more un- 
derstood and tested; in these regions the dependence of the 
line emissivity from the various physical parameters is much 
simpler than in the case of PDRs. As a consequence, lines 
arising mostly from Hii regions are more suitable candi- 
dates to measure the gas metallicity of the inter-stellar mat- 
ter. Therefore, it is very important to understand the emis- 
sivity distribution of various fine-structure emission lines 
in gas clouds and, specifically, understand which ones arise 
mostly from Hii regions and which are instead significantly 
contributed by PDRs. This is possible by using the publicly 
available photoionization model code Cloudy, as described 
by Abel et al. (2005, 2006). 

We have carried out photoionization model calculations 
by using Cloudy version 08.00 (Ferland et al. 1998; Ferland 
2008). We assume pressure-equilibrium gas clouds with a 
plane-parallel geometry. As described later, the results are 
not significantly different if we assume constant-density 
clouds instead of pressure-equilibrium clouds, as far as 
only Hii regions are considered (see §§3.2). Calculations are 
stopped at the depth where the fraction of the H2 molecule 
reaches 50% of the total hydrogen, when we investigate lines 
arising both from Hii regions and from PDRs. This depth 
is defined as the outer edge of PDRs, following Abel et al. 
(2005). On the other hand, calculations are stopped at the 
depth where the fraction of H+ reaches down to 1% of the 
total hydrogen, when the emission lines arising only from 
Hii regions are investigated (§§3.2). We define this depth 
as the outer edge of Hii regions. 

Some low-excitation fine-structure lines can arise be- 
yond this outer edge of PDRs (i.e., fully molecular regions 
surrounding the inner PDRs; see, e.g., Abel et al. 2005). 
We do not investigate such fully-molecular regions since 
our main interest in this work is on emission-lines arising 
within Hii regions, as described below. 

The free parameters in our calculations are (1) the hy- 
drogen density of a cloud (nn) at the irradiated surface of 
a cloud; (2) the ionization parameter (U), i.e., the ratio 
of the ionizing photon density to the hydrogen density at 
the surface of a cloud; (3) the elemental composition of the 
gas; and (4) the spectral energy distribution (SED) of the 
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Table 1. Adopted solar elemental abundance ratios and 
gas-phase depletion factors 
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See the main text for details of the treatments for helium 
and nitrogen. 



photoionizing continuum radiation. We investigate plane- 
parallel gas clouds with gas densities nn — 10^ — 10^ cm~^ 
at the illuminated face (note that the gas density varies 
as a function of the depth within gas clouds), and ioniza- 
tion parameters U = 10~^'^ — 10~^'^. Here the ionization 
parameter is a dimensionless parameter defined as follows: 

u=^-^=^r^d. (1) 

CTlH criH hv 

where c is the speed of the light, ^(H) is the surface flux of 
hydrogen-ionizing photons, F^, is the surface energy flux of 
the input radiation, and vi is the frequency of the Lyman 
edge {hvx = 13.6 eV). 

The adopted solar elemental abundance ratios (in terms 
of number density, not mass) and gas-phase depletion fac- 
tors are given for convenience in Table 1. The adopted de- 
pletion factors arc derived from the literature (e.g., Cowie 
& Songaila 1986; Jenkins 1987; Cardelh et al. 1991; see 
Ferland 2006 for more details). However, there is no broad 
consensus on those factors and, moreover, they may depend 
on the gas density (e.g., Spitzer 1985). We do not discuss 
this issue further in this work, and we simply adopt the de- 
pletion factors given in Table 1. For non-solar metallicities 
we assume that both the dust model and the depletion fac- 
tors are unchanged, but the dust abundance is assiimed to 
scale linearly with the gas metallicity. All elements except 



nitrogen and helium are taken to be primary nucleosynthe- 
sis elements. For helium, we assume a primary nucleosyn- 
thesis component in addition to the primordial value, that 
is, 

— = 0.0737 + 0.0293 Z^^/Zq (2) 
riH 

(Groves et al. 2004). Nitrogen is assumed to be a secondary 
nucleosynthesis element above metallicities of 0.23 solar, 
and a primary nucleosynthesis element at lower metallici- 
ties, i.e., 

log — = 2 log - 4.070 when Z^^ > 0.23^© (3) 
riH Zq 

and 

log — = log ^ - 4.709 when Z^^^ < 0.23Zq (4) 
riH Zq 

(Kewley & Dopita 2002). Note that the behavior of nitrogen 
as a secondary element (Eq. 3) is observationally confirmed 
through the data of Hll regions in galaxies (e.g., van Zee 
et al. 1998; Pilyugin et al. 2003). Chemical evolutionary 
models of galaxies also predict such a behavior for a wide 
range of parameters, including variations of the IMF (e.g., 
Hamann & Ferland 1993; Chiappini et al. 2003). The va- 
lidity of Eq. 3 is sometimes questioned at very high metal- 
licities, Zga,s > 5Zq, which is however beyond the range 
studied in this work (see, e.g., Nagao et al. 2006b; Jiang et 
al. 2008). We investigate gas clouds with gas metallicities 
Zgas = 0.05, 0.2, 0.4, 1.0, 2.0, and 3.0 times solar metallicity. 
In the calculations, Orion-type graphite and silicate grains 
(Baldwin et al. 1991; Ferland 2006), as well as polycyclic 
aromatic hydrocarbon (PAH), are included. 

The adopted ionizing continua are the starburst99 con- 
stant star-formation model spectra (Lcithercr et al. 1999) 
with an age of 1 Myr and an a = 2.35 IMF (Miow = IM© 
and Mup = lOOM©). We selected these model parameters 
because, in this work, we want to simulate actively star- 
forming galaxies (such as infrared galaxies and submillime- 
ter galaxies). However, to demonstrate that our results do 
not depend on the ionizing continuum shape, we will also in- 
vestigate models adopting SEDs with older ages (§3.3). The 
starburst99 spectra are available for stellar metallicities of 
Z^iar = 0.05, 0.2, 0.4, 1.0, and 2.0 times solar metallicity, 
which are shown in Figure 1. For models with Zgg,s = 3.0.^0, 
we adopt the starburst99 spectrum with .^star = ^-OZq. 
The cosmic background emission from radio to X-ray is 
also included (see Ferland 2006 for details). However, this 
component does not affect the calculation results signifi- 
cantly. The heating and ionization by cosmic rays are also 
included, which affects the thermal solutions especially in 
the PDRs (see, e.g., Ferland & Mushotzky 1984; Ferland 
2006). 

3. Results 

3.1. Contributions from H II regions and PDRs 

We first investigate from which parts of clouds emission 
lines originate in order to explore optimal emission-line 
combinations which can be used as good metallicity diag- 
nostics. In Figure 2. we show hydrogen density structures 
and cumulative intensities of some fine-structure emission 
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Fig. 2. Hydrogen density structure (upper panels) and cumulative intensity of some fine-structure emission lines (lower 
panels) of clouds with Zgas = 0.2Zq (left hand panels) and Zgas ~ S.OZq (right hand panels), as a function of hydrogen 
column density. For both cases, nn = 10"^ cm^^ and U = 10^^'^ are adopted (where nn is the hydrogen density at 
the illuminated face of clouds). Solid, dashed, and dotted lines in the upper panels denote the densities of H+, H°, and 
H2, respectively. In the lower panels, solid lines show cumulative intensities for fine-structure emissions mainly from Hll 
regions while dashed lines are for emissions that are emitted also, or mostly from PDRs. 



lines of clouds with Zg^s = 0.2Zq and S.OZq, as func- 
tions of the hydrogen column density. The outer edges of 
PDRs defined above correspond to the end points of calcu- 
lations. The outer edges of Hii regions are clearly seen as 
a sharp drop of the density of H+. It is shown in Figure 
2 that relatively high excitation lines such as [Oiii]51.80 
and [Nlii]57.21 arise within Hii regions while lower excita- 
tion lines such as [Oi]63.17 and [Cii] 157.6 arise mostly in 
PDRs. 

To see how Hii regions and PDRs contribute to the 
intensities of various infrared fine-structure emission lines 
more quantitatively, we give the fraction of the Hll re- 
gion contribution to total emission-line intensity for the 
cases of Zgas = 0.2^0 and 3.0^0, rin = 10^''' cm~^ and 
lO^ " cm--^ and U = lO"^.^ .jnd lO^^ ^ in Tables 2 and 
3. The following fine-structure emission lines arise within 
Hii regions without significant contribution from PDRs (at 
least in the parameter space investigated here): [Neii]12.81, 
[Neiii]15.55, [Siii]18.67, [Siii]33.47, [Neiii]36.01, [Oiii]51.80, 
[Niii]57.21, [Oiii]88.33, [Nii]121.7, and [Nii]205.4 (hereafter 
"Hii region lines"). On the other hand, line intensities of 
the following low-ionization emission lines are significantly 
contributed by the PDRs: [Siii]34.81, [Oi]63.17, [Oi]145.5, 
and [Cii]157.6 (hereafter "PDR lines"). Although the PDR 
lines arise mostly within PDRs, Hii regions also contribute 
to PDR lines in some cases. Note that the contribution by 
Hii regions given in Tables 2 and 3 for those low-ionization 
PDR lines are in some cases upper limits. This is because 
such fine-structure emission lines could be contributed also 



from molecular cloud regions that extend beyond PDRs 
(see, e.g., Abel et al. 2005). 

It is evidently safe to avoid combinations of an Hll re- 
gion line and a PDR line to investigate possible candidate 
of metallicity- diagnostic flux ratios. We therefore focus only 
on Hii region lines to avoid significant parameter depen- 
dences in emission-line flux ratios. 



3.2. Metallicity-sensitive emission- line flux ratios 

Metallicity diagnostics by using infrared fine-structure lines 
have been also discussed in some previous studies. The main 
approach consists in calculating abundances of ionic species 
directly based on their emissivity. The flux (-Fx+> ) of ^ given 
ionic species (X+*) is the product of the density of the ion 
(nx+0) the electron density (rig), and the emissivity of the 
ion (jx+O- Therefore, in Hii regions (where nH+ — nu by 
definition), the following relation applies: 



-Fh+ nY[+ Tic 
which can be re-written as 



riH 



:'H+ 



(5) 



(6) 



where the emissivities are calculated once the gas tem- 
perature is known from the photoionization model. Once 
ionic abundances are inferred by photoionization models. 
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we can estimate the elemental abundances by applying ap- 
propriate ionization correction factors (ICFs). For instance, 
Verma et al. (2003) derived elemental abundances of Ne 
from [Neil] 12.81 and [Neiii] 15.56, Ar from [Arii]6.99 and 
[Ariii]8.99, and S from both [Siii]18.71 and [Siv]10.51, by 
using hydrogen recombination lines such as Bra, for nearby 
starburst galaxies based on Spitzer spectra (see also Wu et 
al. 2008). Panuzzo et al. (2003) also investigated the strat- 
egy to derive N abundance by using [Nii]121.7, [Niii]57.21, 
and Bra. 

However, the requirement of relying hydrogen recombi- 
nation lines makes this approach highly uncertain. Indeed, 
the relatively large wavelength separation between hydro- 
gen recombination lines and metallic emission lines (espe- 
cially for the case of N) is hard to be covered by the same 
instrument, and the use of different instruments or facilities 
may introduce significant systematic errors (due to, e.g., 
aperture effects and calibration errors) in the estimated 
metallicities. Moreover, hydrogen recombination lines are 
sometimes difficult due to the effect of stellar absorption 
features. In addition, the wavelength of Bra is not long 
enough to avoid significant dust extinction in many cases, 
since it is known that the amount of extinction toward the 
most obscured part of star formation in galaxies can reach 
up to Ay ^ 100 mag (e.g., Verma et al. 2003). Although 
Panuzzo et al. (2003) discussed the use of the radio contin- 
uum emission instead of hydrogen recombination lines to 
overcome the above difficulties, the radio continuum gener- 
ally contributed by synchrotron emission, or even an active 
galactic nucleus, and therefore the relation between the ra- 
dio continuum flux and the hydrogen line flux may be com- 
plicated. Therefore, "good" metallicity diagnostics should 
consist of fine-structure emission lines (not of hydrogen re- 
combination lines) at long wavelengths, whose wavelength 
separation is not exceedingly large. 

Liu et al. (2001) investigated N/0 elemental abundance 
ratios of 51 Galactic planetary nebulae (PNs) and proto- 
planetary nebulae (PPNs) with [Oiii]51.80, [Oiii]88.33, and 
[Niii]57.21. They derived the N'^+/0'^+ ratio and trans- 
formed it into the N/0 ratio. Here the two [Oiii] lines are 
used to correct for density effects. The combination of these 
fine-structure emission lines appears to be an optimal set 
of candidate metallicity diagnostics since all of their wave- 
length are long enough, and all of them are simultaneously 
covered by using some observational facilities, such as the 
Herschel Photodetector Array Camera and Spectrometer 
(PACS) (although it requires a small redshift of z > 0.08 
to have all lines in band). We therefore focus on this ap- 
proach and investigate it further below. 

In Figure 3, the predicted [Oiii]51.80/[Niii]57.21 fiux 
ratio is shown as a function of the gas metallicity. It is 
shown that this flux ratio has a strong metallicity depen- 
dence, especially in the high metallicity range (^gas > 
0.2^0). However, this flux ratio also displays (smaller) 
dependences on the ionization parameter and on the 
gas density, which prevent to accurately infer the gas 
metallicity. Figure 4 shows the predicted flux ratios of 
[Oiii]88.33/[Niii]57.21, as a function of the gas metallic- 
ity. Again we wee the density effects on the flux ratio. Note 
that models with a lower gas density predict larger flux ra- 
tios of [Oiii]88.33/[Nlii]57.21, that is in an opposite trend 
from the case of [Oiii]51.80/[Niii]57.21 (Figure 3). 

The diagnostic accuracy improves significantly if the 
sum of [Oiii]51.80 and [Oiii]88.33 is used instead of only 
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Fig. 3. Predicted emission- line flux ratio of 

[Oiii]51.80/[Niii]57.21 as a function of gas metalhcity. 
Blue and red lines denote the models with U = 10~^'^ 
and 10~^-^, and solid and dashed lines denote the models 
with TT-H = 10^'° cm~^ and 10'^'° cm~^, respectively. The 
flux ratio observed in M82 and in the Antennae galaxy is 
shown by black horizontal line (these two galaxies show 
very simflar [Oiii]51.80/[Niii]57.21 flux ratios; see Table 
5). The a;-range of this horizontal line corresponds to the 
inferred metallicity range for M82 and the Antenna galaxy. 
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Fig. 4. Same as Figure 3 but for the flux ratio of 
[Oiii]88.33/[Niii]57.21. 



LO 



10 r 



00 

? 1 

o 

00 



o 



— 1 — I — I — [—r 




logn=1 .0 
logn=3.0 



logi7=-1.5 logC/=-2.5 



0.1 1 
^qas / ■^O 



Fig. 5. Same as Figure 3 but for the fiux ratio of 
([Oiii]51.80-l-[Oiii]88.33)/[Niii]57.21. Note the much lower 
dependence on the gas density, making this ratio particu- 
larly suited to measure the gas metallicity. 



6 



Nagao et al.: Infrared Metallicity Diagnostics 



in 1 Q 
in ' ^ 

iri 



CD 



00 

^ 0.1 

z 

0.01 



- logn=1 .0 



(N 

in 



Csl 



0.1 



0.01 



r-r-r-^ 

logn=2.0 



logn=3.0 



"n 

logC/=-3.5- 




logC/=-2.0" 
r logi7=-1.& 

^\ogU=-2.5'- 

:\ogU=-3.0 
• logC/=-2.0 

'- \ogU=-^.5 



-. logC/"=-2.0 



logn=3.0 



logC/=-1.5- 



logn=4.0 




0.5 1 2 5 10 

[0III]51.88/[0III]88. 33 

Fig. 6. Predicted emission-line flux ratio of 
[Neii]12.81/[Neiii]15.55 (upper) and [Nii]121.7/[Niii]57.21 
(lower) as a function of the density-sensitive flux ratio 
[Oiii]51.80/[Oiii]88.33. Solid and dotted lines denote 
[/-constant and rt-constant model sequences, respectively. 
Models with Zg^s = 0.2, 1.0, and 3.0 Zq are denoted by 
blue, green, and red lines, respectively. The models with 
Zgas = 0.2Zq or with rin = 10^ ° cm~^ are not shown 
in the lower panel because [Nil] 121.7 is extremely faint 
in those case. Observational data for M82 and for the 
Antennae galaxy are shown in the lower panel. 



[Oiii]51.80, as shown in Figure 5. The density effect on the 
flux ratio of ([Oiii]51.80+[Oiii]88.33)/[Niii]57.21 is smah 
(^ 50%) over the whole range of metallicities. However, 
there is still a significant dependence on the ionization 
parameter, especially at high-metallicities {Zg^s ^ 2Zq), 
where it may cause variation of the ratio up to a factor of 
~ 3. We will compare these model predictions with obser- 
vational data in §§4.1. 

The predicted fluxes of Hii region lines are summarized 
in Table 4, where they are normalized by the [Oiii]51.88 
flux. The results do not change significantly when we adopt 
constant-density models instead of constant-pressure mod- 
els. More specifically, the differences in the predicted fluxes 



of the emission lines given in Table 4 between constant- 
pressure models and constant-density models are less than 
20% in most cases. 

Next we focus on how to account for the dependence of 
the ionization parameter of the flux ratios shown in Figures 
3, 4, and 5. To estimate the ionization parameter, pairs of 
ion species of the same element in a different ionization 
degree have been often used. We thus examine two flux ra- 
tios of [Neii]12.81/[Neiii]15.55 and [Nii]121.7/[Niii]57.21 as 
possible tools to correct for the ionization parameter effect. 
Note that also [Neii]12.81 and [Nii]121.7 are Hii region hues 
(i.e. not from PDRs), as described in §§3.1. In Figure 6, 
we show these two emission-line flux ratios as a function 
of the density-sensitive flux ratio [Oiil]51.80/[Oiii]88.33. 
Clearly, both of the two flux ratios are strongly depen- 
dent on the ionization parameter and therefore they can 
be used to measure U and correct the metallicity diagnos- 
tic flux ratio ([Oiii]51.80+[Oiii]88.33)/[Niii]57.21 for this 
dependence. These two flux ratios are also sensitive to 
metallicity, that is partly due to the metallicity-dependent 
SED of the ionizing photons (Figure 1). Therefore an iter- 
ative process is required to estimate metallicity of galax- 
ies by using the diagnostic diagrams given in Figures 5 
and 6. Which is the [/-sensitive flux ratio more appro- 
priate to correct for the ionization parameter effects on 
the metallicity determination depends on the speciflc case. 
The [Neil] 12.81 and [Neiii] 15.55 emission lines are mod- 
erately strong in star-forming galaxies and have been ob- 
served by Spitzer in many galaxies (e.g.. Dale et al. 2006), 
and their flux ratio has only a slight dependence on the 
gas density, as shown in Figure 6. However, these lines 
are more affected by dust extinction effect with respect 
to the metallicity diagnostic emission lines, [Oiii]51.80, 
[Oiii]88.33, and [Niii]57.21. In contrast, [Niii]57.21 and 
[Nii]121.7 do not suffer significant dust extinction. However, 
the expected emissivity of [Nii]121.7 is low, which makes 
its measurement challenging. In Figure 7, we also plot 
the same [/-sensitive flux ratios, [Neii]12.81/[Nelll]15.55 
and [Nii]121.7/[Niii]57.21 as a function of the density- 
sensitive flux ratio [Nii]121.7/[Nii]205.4, instead of 
[Oiii]51.80/[Oiii]88.33. Since the [Nii]205.4 emission is very 
faint and at a very long wavelength, it is very difficult 
to study with Herschel and SPICA, though its detection 
should be feasible with ALMA. 



3.3. Dependences on the stellar age 

All emission-line flux ratios shown in Figures 3-7 are cal- 
culated by adopting constant star-formation SEDs with an 
age of 1 Myr. However, it should be verified whether the 
predicted fiux ratios depend on the age of the stellar popu- 
lation since the adopted age (1 Myr) seems too young com- 
pared to the typical age of star-forming galaxies in general. 
Cid Fernandes et al. (2003) investigated stellar populations 
of nearby starburst galaxies and reported that the typical 
starburst age is ~ iqT-O-t.b yj._ More recently, Rodriguez 
Zaurin et al. (2010) studied stellar populations of low-z 
ULIRGs, finding similar age ranges. It is thus important to 
examine how the emission-line diagnostics studied in this 
work depend on the age of the stellar populations used for 
the input SED. 

We calculated the emission-line flux ra- 
tios of ([Oiii]51.80-h[Oiii]88.33)/[Niii]57.21 and 
[Nii]121.7/[Niii]57.21 by adopting the input SED of 
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Fig. 7. Predicted emission-line flux ratio of 
[Neii]12.81/[Neiii]15.55 (upper) and [Nii]121.7/[Niii]57.21 
(lower) as a function of the density-sensitive flux ratio 
[Nii]121.7/[Nii]205.4. Solid and dotted lines denote U- 
constant and nu-constant model sequences, respectively. 
Models with Zgas = 1-0 Zq and 3.0 Zq are denoted by 
green and red lines, respectively. Only the models with 
riH < 10^ cm~^ are plotted for the cases of Zg^s = I.OZq, 
and models with nn < 10'^ cm~^ are plotted for the 
cases of Zgas = S.OZq. This is because the [Nii] lines are 
faint in clouds with a high density due to the collisional 
de-excitation effects. 



stellar populations with ages of 3, 10, and 30 Myr in 
addition to our default value (1 Myr). Figure 8 shows the 
results of the age dependences of those two diagnostics. 
The flux ratio of ([Oiii]51.80-H[Oiii]88.33)/[Niii]57.21 
changes by less than 30% in the range 1 — 30 Myr, that 
is negligible with respect to the metallicity dependence of 
this flux ratio shown in Figure 5. The age dependence is 
not significant also for the flux ratio [Nii]121.7/[Niii]57.21. 

These results are due to the fact that the SED shape 
significantly changes as a function of the stellar-population 
age (for the range of 1 — 30 Myr) only at A > 3000A (see 
Figure 8 in Leitherer et al. 1999), when the constant star- 
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Fig. 8. Dependence of the predicted emission-line fiux ra- 
tios ([Oiii]51.80+[Oiii]88.33)/[Niil]57.21 (upper panel) and 
[Nii]121.7/[Niii]57.21 (lower panel) on the age of the stel- 
lar population used in the photoionization model calcula- 
tions. Models with constant star-formation history and with 
solar metallicity are displayed. The meanings of lines are 
the same as in Figure 3. In the lower panel, models with 
TiH = 10'^ cm~^ and U — 10^^'^ are not plotted, because 
the flux of the [Nii]121.7 is too faint for those parameters. 
The displayed ranges of emission-line flux ratios in the up- 
per and lower panels are the same as those in Figures 5 and 
the lower panel of Figure 6, respectively. 



formation history is assumed. The SED at shorter wave- 
lengths does not significantly depend on the starburst age, 
because in that wavelength range it is dominated by emis- 
sion from massive stars. Since the ionization structure in Hii 
regions is mostly determined by ultraviolet photons (whose 
SED is not sensitive to the stellar-population age) , the flux 
ratios of Hii-region lines are consequently insensitive to the 
starburst age. These situation may be completely different 
for FDR lines such as [Cii] 157.6, but the age dependence of 
such PDR lines is beyond the scope of this work. Note that 
the constant star-formation history is a better assumption 
than the single-burst model for star-forming galaxies when 
relatively old ages (^ 30 Mye) are assumed, because stellar 
populations with an age of ~ 30 Myr in the single-burst 
model do not produce ionizing photons and thus we do not 
expect detectable emission-line fluxes from such systems. 
It can be therefore concluded that our analysis on metal- 
licity (and other) diagnostics does not significantly depend 
on the age of stellar populations. 
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4. Discussion 

4.1. Comparison with observational data 

Here we compare our calculation results on long- wavelength 
(A > 50/im) diagnostic fine structure lines with existing ob- 
servational data. The only previous instrument that mea- 
sured the diagnostic fine-structure lines in star-forming 
galaxies is LWS onboard ISO. Specifically we focus on 
two galaxies, M82 and the Antennae galaxy, since many 
infrared fine-structure lines are measured for them (Fischer 
et al. 1996; Colbert et al. 1999). The measured fluxes 
are summarized in Table 5. Interestingly, the measured 
flux ratio of [Oiii]51.80/[Oiii]88.33 is quite similar be- 
tween those two galaxies; 1.2 and 1.0 for M82 and the 
Antennae. These fiux ratios suggest nn ~ 10^ cm~^ for U ^ 
10~^ (Figure 5). The metallicity diagnostic fiux ratio of 
([Olil]51.80-H[Oiii]88.33)/[Niii]57.21 is also similar between 
M82 and the Antennae, 5.6 and 6.0 respectively. Taking the 
inferred gas density also into account, these fiux ratios cor- 
respond to Z ^ 0.5 — 0.7.^0, i.e., slightly sub-solar metal- 
licities. The observed flux ratio of [Nii]121.7/[Niii]57.21 in 
M82 is measured to be 0.5, which, combined with the our 
initial metallicity estimation, gives logU = —3.5. Only an 
upper limit of [Nii]121.7/[Niii]57.21< 0.28 (3 a) is reported 
for the Antennae galaxy, providing an upper limit on log U 
of about —3. Although the constraint on the ionization pa- 
rameter is not strong, the uncertainty on the ionization 
parameter does not affect significantly the metallicity esti- 
mation. 

Although both M82 and the Antennae are well-studied 
nearby galaxies, it is not so straightforward to compare 
the inferred gas metallicity based on the fine-structure 
emission-line flux ratios with that inferred from other 
metallicity diagnostics, because the aperture size of LWS 
(the beam-size is ~ 80 arcsec in FWHM) is very large and 
difficult to compare with other measurements. Origlia et al. 
(2004) measured the nuclear gaseous and stellar metallici- 
ties by using X-ray and near-infrared spectroscopic observa- 
tions, and they found that both the metallicities are close to 
or slightly less than the solar metallicity in M82. Taking the 
metallicity gradient into account, these metallicities (based 
also on diagnostics little affected by dust extinction) are 
consistent with the metallicity inferred through the far-IR 
fine-structure diagnostics. The situations is similar also for 
the Antennae galaxy, since its gaseous and stellar metallic- 
ities are close to or slightly less than the solar metallicity 
(Bastian et al. 2006, 2009). Although the comparisons dis- 
cussed above do not provide a tight test for our new metal- 
licity diagnostics based on far-IR fine structure lines, do 
at least suggest a broad consistency with other metallicity 
diagnostics. 

4.2. Observational feasibility with next-generation 
instruments 

We finally examine the detectability of the diagnostic 
fine-structure emission lines discussed here with the new- 
generation instruments, PACS on Herschel (Poglitsch et 
al. 2010) and SAFARI (Swinyard 2008; Swinyard et al. 
2009) the spectrometer planned for SPICA (a 3m-class 
cooled telescope in space, Nakagawa 2009). Figure 9 shows 
the detectability of [Nii]121.7, [Niii]57.21, [Oiii]51.80, and 
[Oiii]88.33 with PACS. Specifically the figure shows the 
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Fig. 9. Detectability with PACS on Herschel of some di- 
agnostic far-IR fine-structure lines discussed in this paper, 
which can be used to measure the gas metallicity, ionization 
parameter and density in galaxies, as a function of redshift. 
The left hand scale gives minimum LpiR to achieve a 5a de- 
tection of each line with an exposure time of 3 hours. The 
right hand axis shows the minimum star formation rate for 
the line detectability corresponding to the minimum LpiR, 
based on the relation given in Kennicutt (1999). Green, 
red, blue, and magenta lines denote [Nii]121.7, [Niii]57.21, 
[Oiii]51.80, and [Oiii]88.33, respectively. 
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Fig. 10. Same as Figure 9 but for the detectability with 
SAFARI boarded on SPICA. 
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mininium far-infrarcid luminosity of galaxies to have each 
far-IR line detected (S/N=5) with a 3 hour exposure, 
as a function of redshift. Here we make the simplify- 
ing assumption that the luminosity of these emission 
lines scales roughly with far-infrared luminosity and keep- 
ing a ratio similar, on average, to what observed in 
some nearby galaxies where all of these transitions were 
observed with past ISO observations; more specifically, 

-^[OIII]88.33/-^FIR ~ 1-8 X 10~^, i[oiII]51.8o/^FIR ~ 1-8 X 
10"^, -£'[NIII]57.2l/-f'FIR ~ 6.0 X 10"'', and I'[NII]121.7/-f'FIR ~ 

2.7 X 10"^ (Fischer et al. 1996; Colbert et al. 1999). The 
scale on the right hand side; gives the minimum SFR cor- 
responding the minimum LpiR, by using the SFR-Lfir 
relationship given in Kennicutt (1999). Figure 9 suggests 
that we can detect all the fine-structure emission lines 
needed to infer the metallicity, in luminous infrared galaxies 
(Lfir ^ 10" - lO^^L©) at rcdshifts up to z - 0.4. The de- 
tection to somewhat higher redshift is restricted to ULIRGs 
(Lfir > 10^^ Lq). On the other hand, SAFARI boarded on 
SPICA will drastically expand the detectability, as shown in 
Figure 10. By using SAFARI, the diagnostic fine-structure 
lines will be detectable for star-forming galaxies at z ~ 1 or 
at even higher redshift. In Figure 10 we adopt the SAFARI 
sensitivity given by Swinyard (2008) with a scaling of sen- 
sitivity for a conservative 3.0m aperture. Actually the red- 
shift limit may be limited not by the instrument sensitivity, 
but by the wavelength coverage, although the instrument 
design of SAFARI has not yet completely fixed. 

5. Summary 

Gas metallicity diagnostics with infrared fine-structure 

emission lines arc investigated in this paper, based on de- 
tail photoionization model calculations. The main results 
are as follows: 

— Some fine-structure lines arise both in Hll regions and 
PDRs and thus only emission lines arising within Hll re- 
gions (not in PDRs) should be used as good metallicity 
diagnostics. 

— Emission-line flux ratios that consist of [Oiii]51.80, 
[Nlli]57.21, and [Olll]88.33 provide good metal- 
licity diagnostics; especially the flux ratio of 
([Oiii]51.80-F[Oiii]88.33)/[Niii]57.21 is sensitive to 
the gas metallicity with a small dependence on the gas 
density. 

— To correct for the effects of the ionization parame- 
ter on the infrared metallicity diagnostics, flux ratios 
of [Ncii]12.81/[Nciii]15.55 and [Nii]121.7/[Niii]57.21 are 
useful. 

— The detection of those infrared fine-structure lines is 
feasible, up to z ^ 0.4 with using Herschel/PACS and 
even up to z ^ 1 or more with using SPIC A/SAFARI. 

Since the infrared fine-structure emission lines are not af- 
fected by the dust extinction, the gas-mctallicity diagnos- 
tics presented in this paper are usehil to assess the chemical 
properties in dust-obscured populations such as infrared 
galaxies, that are very important in terms of the galaxy 
evolution. Some of the infrared fine-structure emission lines 
from high-z galaxies will be apparently exciting targets also 
for ALMA. 
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Table 2. Fractions of the Hii region contribution to total emission- line fluxes in low-metallicity cases {Zga,s = 0-2Zq) 
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Table 3. Fractions of the Hii region contribution to total emission- line fluxes in high-metallicity cases (^gag — 3Zq^ 



Line 


riH = 10^ ° cm"* 

u = 10-2-^ 


TiH = 10^ ° cm"^ 

u = 10-^-^ 


riH = 10^ ° cm"^ 

u = 10-2-^ 


riH = 10^ ° 

u = 10-1-^ 


[Ne ii] 12.81 


1.00 


1.00 


1.00 


1.00 


[Ne III] 15.55 


1.00 


1.00 


1.00 


1.00 


[S III] 18.67 


1.00 


1.00 


1.00 


1.00 


[S III] 33.47 


1.00 


1.00 


1.00 


1.00 


[Si II] 34.81 


1.00 


0.96 


0.61 


0.37 


[Nc III] 36.01 


1.00 


1.00 


1.00 


1.00 


[0 III] 51.80 


1.00 


1.00 


1.00 


1.00 


[N III] 57.21 


1.00 


1.00 


1.00 


1.00 


[0 I] 63.17 


0.57 


0.06 


0.03 


0.02 


[0 III] 88.33 


1.00 


1.00 


1.00 


1.00 


[N ii] 121.7 


1.00 


1.00 


1.00 


1.00 


[O i] 145.5 


0.82 


0.12 


0.05 


0.02 


[C II] 157.6 


0.75 


0.34 


0.15 


0.11 


[N ii] 205.4 


1.00 


1.00 


1.00 


1.00 
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Table 4. Predicted emission-line fluxes by photoionization models 



Model Parameter Predicted Fluxes Normalized by F([Oiii]51.88) 

log nu'' log U Zge^/ZQ [Neii]12.81 [Neiii]15.55 [Oiii]51.88 [Nni]57.21 [Oiii]88.33 [Nii]121.7 [Nii]205.4 



1.0 


-1.5 


0.05 


0.004 


0.809 


1.000 


0.106 


1.664 


<0.001 


<0.001 


1.0 


-1.5 


0.20 


0.006 


0.789 


1.000 


0.145 


1.646 


<0.001 


<0.001 


1.0 


-1.5 


0.40 


0.011 


0.769 


1.000 


0.274 


1.617 


0.002 


0.002 


1.0 


-1.5 


1.00 


0.054 


0.699 


1.000 


0.698 


1.536 


0.009 


0.008 


1.0 


-1.5 


2.00 


0.394 


0.396 


1.000 


1.845 


1.488 


0.078 


0.055 


1.0 


-1.5 


3.00 


0.618 


0.287 


1.000 


3.147 


1.471 


0.314 


0.181 


1.0 


-2.5 


0.05 


0.031 


0.928 


1.000 


0.171 


1.691 


0.006 


0.006 


1.0 


-2.5 


0.20 


0.053 


0.916 


1.000 


0.179 


1.693 


0.007 


0.007 


1.0 


-2.5 


0.40 


0.089 


0.899 


1.000 


0.320 


1.698 


0.014 


0.015 


1.0 


-2.5 


1.00 


0.322 


0.778 


1.000 


0.937 


1.708 


0.073 


0.074 


1.0 


-2.5 


2.00 


1.382 


0.409 


1.000 


3.229 


1.737 


0.598 


0.543 


1.0 


-2.5 


3.00 


1.916 


0.328 


1.000 


5.308 


1.783 


1.714 


1.423 


3.0 


-1.5 


0.05 


0.004 


0.942 


1.000 


0.042 


0.408 


<0.001 


<0.001 


3.0 


-1.5 


0.20 


0.008 


0.984 


1.000 


0.058 


0.362 


<0.001 


<0.001 


3.0 


-1.5 


0.40 


0.015 


1.039 


1.000 


0.109 


0.329 


<0.001 


<0.001 


3.0 


-1.5 


1.00 


0.080 


1.132 


1.000 


0.275 


0.284 


<0.001 


<0.001 


3.0 


-1.5 


2.00 


0.529 


0.825 


1.000 


0.612 


0.257 


0.006 


<0.001 


3.0 


-1.5 


3.00 


0.767 


0.787 


1.000 


0.924 


0.243 


0.015 


0.002 


3.0 


-2.5 


0.05 


0.036 


1.057 


1.000 


0.068 


0.417 


<0.001 


<0.001 


3.0 


-2.5 


0.20 


0.062 


1.061 


1.000 


0.071 


0.401 


<0.001 


<0.001 


3.0 


-2.5 


0.40 


0.106 


1.067 


1.000 


0.127 


0.383 


0.003 


<0.001 


3.0 


-2.5 


1.00 


0.395 


0.986 


1.000 


0.364 


0.348 


0.013 


0.002 


3.0 


-2.5 


2.00 


1.607 


0.553 


1.000 


1.094 


0.316 


0.069 


0.009 


3.0 


-2.5 


3.00 


2.245 


0.534 


1.000 


1.757 


0.298 


0.155 


0.019 



Total hydrogen density at the illuminated surface, in units of cm 



Table 5. Measured fluxes of infrared fine-structure lines in M82 and the Antennae galaxy* 



line 


M82'' 




Antenna" 


[Oiii]51.80 


10.3 ± 0.5 X 10' 


-18 


4.9 ±0.29 X 10"^^ 


[Niii]57.21 


3.4 ±0.5 X 10" 


-18 


1.6 ±0.058 X 10"^" 


[Oiii] 88.33 


8.6 ± 0.4 X 10" 


-18 


4.7 ±0.13 X 10"^^ 


[Nil] 121.7 


1.7 ±0.3 X 10" 


-18 


< 0.44 X 10-1^ (3ct) 



Line fluxes arc given in units of W cm 
Data taken from Colbert et al. (1999). 
" Data taken from Fischer et al. (1996). 



